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Abstract 

We investigate the hypothesized existence of an S-matrix for gravity, and some of its 
expected general properties. We first discuss basic questions regarding existence of such 
a matrix, including those of infrared divergences and description of asymptotic states. 
Distinct scattering behavior occurs in the Born, eikonal, and strong gravity regimes, and 
we describe aspects of both the partial wave and momentum space amplitudes, and their 
analytic properties, from these regimes. Classically the strong gravity region would be 
dominated by formation of black holes, and we assume its unitary quantum dynamics is 
described by corresponding resonances. Masslessness limits some powerful methods and re- 
sults that apply to massive theories, though a continuation path implying crossing symme- 
try plausibly still exists. Physical properties of gravity suggest nonpolynomial amplitudes, 
although crossing and causality constrain (with modest assumptions) this nonpolynomial 
behavior, particularly requiring a polynomial bound in complex s at fixed physical momen- 
tum transfer. We explore the hypothesis that such behavior corresponds to a nonlocality 
intrinsic to gravity, but consistent with unitarity, analyticity, crossing, and causality. 
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1. Introduction 

In a quantum-gravitational theory where spacetime, locality, etc. may not be funda- 
mental concepts, an important question is what quantities are amenable to quantitative 
analysis. In this paper, we will assume that flat space, or something which it closely ap- 
proximates, is an allowed configuration of the theory. We will moreover assume that there 
is an action of its symmetry group, namely the Poincare group, both on this configura- 
tion and on perturbations about it. This suggests that we can consider states incident 
from infinity, with given momenta and energies, and study their scattering. The resulting 
quantum amplitudes should be summarized in an S-matrix. 

One would like to understand what properties are expected of such an S-matrix. For 
a quantum theory, unitarity is a given. Analyticity in momenta and crossing symme- 
try encode important physical features of S-matrices in quantum field theory (QFT), like 
causality Gravity offers some new features whose role needs to be understood. Mass- 
lessness is first, and causes infrared singularities; these we however envision regulating by 
working in spacetime dimension D > 4, or by proper formulation of inclusive amplitudes. 
Another is growth of the range of gravity with energy, as is seen for example in growth of 
the Schwarzschild radius of a black hole formed in a high-energy collision. An important 
question is how these new features can be reconciled with the others. One would also like to 
understand how these and other physical properties either do or don't manifest themselves 
in a gravitational S-matrix - particularly locality and causality. The latter properties are 
especially interesting, given that a certain lack of locality could be part of a mechanism 
for information to escape a black hole, and thus explain the mysteries surrounding the 
information paradox. Yet locality is manifest in low-energy descriptions of nature, and is 
a cornerstone of QFT; it is also nontrivially related to causality, which plays an important 
role in consistency of a theory. 

In this paper, we carry out some preliminary investigation of these matters, with 
particular focus on the ultra-high energy regime. We will make the maximal analyticity 
hypothesis where one assumes that the only singularities that appear in the scattering 
amplitudes are those dictated by unitarity. Our investigations will then focus on the 
question of what can be learned by combining unitarity, analyticity, crossing and causality 
together with expected general features of gravity. In spite of the plausibly nonlocal 
behavior of the gravitational amplitudes that we will explore, we have found no evidence 
for a lack of harmony between such nonlocality and these basic properties. We thus 
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entertain the possibility that an S-matrix representation of such nonlocal dynamics exists, 
which retains the essential physical features. 

The next section will further describe the S-matrix hypothesis, and some issues that 
must be confronted in its formulation, particularly questions of infrared divergences and 
asymptotic completeness, and summarizes aspects of exclusive amplitudes and their partial 
wave expansion. Section three contains a summary of the different scattering regimes 
(broadly, Born, eikonal, and strong gravity), and aspects of the physics of each. Section 
four focusses on the strong gravity regime, where one expects significant contributions from 
processes classically described as black hole formation. We parameterize the corresponding 
intermediate states as resonances, and investigate their implications for the form of the 
partial wave amplitudes. Section five further develops the description of these amplitudes, 
summarizing our knowledge of the contributions to the phase shifts and their imaginary 
parts from the different regimes. Section six overviews some properties of amplitudes in 
momentum space, some of which can be inferred from those of partial wave amplitudes. 
In particular, for both forms of amplitudes, we find strong indications of non-polynomial 
behavior. Section seven investigates aspects of analyticity and crossing; the latter is less 
transparent than in a theory with a mass gap. Nonetheless, there is an argument for 
crossing, and this together with causality (plus hermitian analyticity and a smoothness 
assumption) in turn leads to constraints on non-polynomial growth. Section eight closes 
with further discussion of nonpolynomiality, and its connection with the question of locality 
of the theory. 

Study of ultraplanckian collisions in gravity has a long history. In string theory, this 
includes [|2|-[5| and |§, and other prominent early references are An important 

question is whether string theory resolves the puzzles of this regime. In particular, the 
information paradox suggests a breakdown of locality in this context; while string theory 



is apparently nonlocal due to string extendedness, it has been argued pllJl^l that this effect 
does not appear to enter in a central way in the regimes of interest. In fact, the strong 
gravitational regime, where classically black holes form, apparently corresponds to a break- 



down of the gravitational loop expansion. Ref. []T3[ has argued for a possible resummation 
of string amplitudes that continues into this regime, but we view the apparent need for 
nonlocal mechanics as well as the absence of clearly relevant stringy effects as suggesting 
that a new ingredient is instead required for fundamental description of this regime |TT|j . 
Though a perturbative string description appears insufficient for a complete description, 
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it has been argued that non-perturbative dual formulations such as AdS/CFTjjnj] will ad- 
dress these problems. While there has been some progress towards extracting a flat space 
gravitational S-matrix from AdS/CFT |]T5| - |T9[| , some puzzles remain [[PTl i pOj about whether 
this is possible; one expects similar issues in Matrix theory plfl . Whether or not it is, we 
take a more general viewpoint, extending work of |[22|| : whatever theory provides this S- 
matrix, we would like to characterize its features, and some of those may be rather special 
in order to describe gravity. Moreover, it may be that, as suggested in P^| , the need to 
describe such features is in fact a critical clue to the dynamics of a quantum theory of 
gravity. 



2. The hypothesis of the gravitational S-matrix 

It is natural to expect that the problem of high-energy gravitational scattering in 
asymptotically flat space can be properly formulated in terms of the S-matrix. Here, 
however, one must grapple with some preliminary issues. 

A first issue is that we don't know a precise description of the quantum numbers of 
these states. For example, they could be states of string theory, some other completion 
of supergravity, or some other theory of gravity. However, in any case, we expect that 
the asymptotic states include those corresponding to widely separated individual incident 
particles, e.g. electrons, neutrinos, etc., in order to match our familiar description of 
nature. Or, if the theory were string theory, incident states are string states. We might have 
states with other quantum numbers as well. An example of the latter that is sometimes 
useful to consider is scattering in Minkowski space that is reached by compactification from 
higher-dimensions; there, one may have incident particles or strings with conserved Kaluza- 
Klein charge. In any of these cases, a nice feature of gravity is that it universally couples 
to all energy, so we view it as plausible that some important features of gravitational 
scattering, particularly at high-energy, are independent of this detailed description of the 
asymptotic states. 

A second issue is that, in a perturbative description of gravitons propagating in flat 
space, gravity suffers from infrared divergences in four dimensions, arising from soft gravi- 
tons, and as a consequence one must generalize from the S-matrix to inclusive amplitudes. 
While it does not seem inconceivable that this is of fundamental importance, we will as- 
sume that it is not. One reason for this is that QED suffers a similar problem, with 
the simple resolution through inclusive generalization of the S-matrix, summing over soft 
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photon states. Moreover, we note that this problem is not present if one works with higher- 
dimensional gravity. Specifically, for spacetime dimension D > 4, soft graviton divergences 
are not present. (For D > 7, the total cross-section is finite.) We have already motivated 
considering higher-dimensional theories, by including the possibilities of string theory or 
supergravity, or we may simply think of this as dimensional regularization - in any case, 
to avoid this issue we will typically work in D > 4. 

Another issue that plausibly comes closer to being fundamental regards the question 
of asymptotic completeness of states. The asymptotic completeness condition^ states that 
the Hilbert space of the theory is equivalent to a Fock space of asymptotic free particles. 
However, there are apparent limitations to such a Fock space description. An example 
is the locality bound []25|j2^ , p6| and its N-particle generalizations |TT| . Specifically, if one 
considers two particles in wavepackets, which we for example can take to be gaussian with 
central positions and momenta x, y and p, q, these have a field theory description in terms 
of a Fock space state <j>x,p<j> y ,q\0) ■ However, such a description must break down when we 
violate the bound 

\x-y\ D ~ 3 > G\p + q\, (2.1) 

where G ~ Go, the D-dimensional Newton constant. In this regime, gravity becomes 
strong, and so limits a Fock space description of the system; this limitation in principle 
extends to arbitrarily large distances. One may yet be able to construct an asymptotic 
description of all states in terms of free-particle states, using further evolution - if one 



evolves a state violating (|2.1|) backwards in time, it generically ceases violating the bound, 
and would be expected to resolve itself into well-separated free particles asymptotic from 
infinity. Thus, with such a limiting procedure, and a weak form of local Lorentz invariance 
(in order to describe asymptotic particles with relative boosts), one plausibly describes 
asymptotics in terms of certain Fock space states. 

In short, we will hypothesize the existence of a gravitational S-matrix, or its inclusive 
generalization in D = 4. While we do not have a complete description of the asymptotic 
states, we will assume that they include states closely approximating particles that are 
initially widely-separated, and moreover are allowed to have very large relative momenta. 
This starting point amounts to making certain assumptions about a weak notion of locality 
(asymptotically separated particles) and local Lorentz invariance (large relative boosts 



1 See, e.g., chapter 7 of 



allowed for widely separated particles). However, we will not necessarily assume that 
stronger forms of locality and local Lorentz invariance are fundamental in the theory. 

For practical purposes, it is often convenient to imagine that the asymptotic states 
correspond to spinless particles of mass m, plus gravitons. With such a collection of 
asymptotic states \a)i n , \ot) ou t, (taken to be Heisenberg-picture states) we expect an S- 
matrix of the form 

S a/3 = ou t(a\/3)in = (a\S\f3) . (2.2) 
As usual, we separate off the non-trivial part as S = 1 + VT . 

2.1. Exclusive amplitudes 

Much of this paper's discussion will focus on the simplest non-trivial amplitude of the 
theory, that for exclusive 2 — ^ 2 scattering. Here, the transition matrix element T (in the 
plane wave limit) is then defined by 

(P3,P4|T|pi,p 2 ) = % 3P a,piv2 = ( 27 D £ D (Pi +P2 -Pa -p4)T(s,t) , (2.3) 

and is a function of the Mandelstam parameters 

8 = -(^ + p 2 ) 2 = E 2 , t = -(pi - p 3 ) 2 , u = -(p! - p 4 ) 2 . (2.4) 

We expect that important features of the theory are encoded in this amplitude and 
its analyticity properties. Since the graviton is massless, amplitudes are singular at t = 0, 
and likewise in other channels; for example, the Born approximation to t-channel exchange 
gives 

T tree (s,t) = -8ivG D s 2 /t. (2.5) 
We will consider other aspects of analyticity in section seven. 

2.2. Partial wave expansion 

Unitarity and some other physical features of the amplitude are most clearly formu- 
lated by working with the D-dimensional partial wave expansion, which is0 

oo 

T(M) = ^ 2 " D/2 ]T(Z + A)QVs#)/z(s) • (2.6) 

1=0 

Here A = (D - 3)/2, 

V* = 2 4A+3 7r A r(A) , (2.7) 
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and are Gegenbauer polynomials, with arguments given by the center-of-mass (CM) 
scattering angle, 

2t 

cos# = l + — r . (2.8) 

s - Am 2 

Note that 

t = (Am 2 - s) sin 2 (#/2) , u = (Am 2 - s) cos 2 (#/2) . (2.9) 



The inverse relationship to ( |2.6| ) gives the partial wave coefficients fi(s) in terms of the 
matrix element, 

(D-4)/2 j-K 

fi(s)= „ Am / d6 S m D - 3 9Cf(cos9)T[s,(Am 2 -s)sm 2 (9/2)} , (2.10) 

iD^i (1) Jo 

with 

7D = 2r(^)(167r)^- 2 )/ 2 . (2.11) 

The unitarity condition 

Imfi(s) > {Ms)] 2 , (2.12) 

for real s > can be solved in terms of two real parameters, the phase shift Si(s), and the 
absorptive coefficients (3i(s) > 0: 

/,( a ) = 1 h_ e 2^( S )-2/3 ! ( S )l . (2.13) 

It is important to understand the convergence properties of the partial wave expansion 
( |2.6|) . For a theory with a mass gap, the expansion can be shown to converge in the 
Lehmann ellipse|27], which extends into the unphysical regime t > 0, cos# > 1. This 
extension is useful for further constraining amplitudes, e.g. through the Froissart-Martin 
H,|9[ bound. 

Masslessness of gravity alters this behavior. Let us first ask when the partial wave 
coefficients (|2.10|) are well defined. Specifically, at long-distance/small angle, we have the 
Born approximation, ( |2.5|) . This gives a pole at zero angle, T ~ I/O 2 , and correspondingly 
the integral ( |2.1U| ) only converges for D > A. While other long-distance effects, like soft 
graviton emission, could modify the amplitude (P75|), we don't expect them to alter this 
convergence behavior. 

In general, a series of the form ( |2.6| ) converges in an ellipse with foci at cos# = ±1. 
The existence of the singularity in T at 9 = indicates that the partial wave expansion 
does not converge past cos# = 1. Thus, the Lehmann ellipse has collapsed into a line 
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segment along the real axis. Note that one does expect ImT(# = 0) to be finite for D > 7. 
This follows from the optical theorem (see the Appendix) - as we have noted, the Born 
cross section given by ([275|) is not infrared divergent for D > 7. However, this finiteness 
does not indicate that the expansion of ImT can be continued past this point - higher 
derivatives of ImT are expected to in general diverge at 9 = 0. 

The failure of convergence of the partial wave expansion in the regime t > is an 
impediment to using some of the powerful methods that have been successfully applied in 
theories with a mass gap. Nonetheless, we suggest that study of partial wave amplitudes 
can still be useful for inferring features of scattering. While we are in particular interested 
in features of the analytic continuation of T(s,t) to complex values of s and t, where 
convergence of the expansion is problematic, we can exploit the inverse relation ( |2.10|) . 
Regardless of the convergence of the partial wave expansion, we have argued that ( |2.10| ) 
is convergent for D > 5. Thus, if physical considerations imply statements about the 
behavior of fi(s), these in turn imply properties of the integrand of ( |2.10| ), and specifically 
of T(s,t). 

3. Scattering regimes 

In different regions of s and t, or E and /, we expect differing physical behavior of 
amplitudes. A more pictorial way to think of these different regimes is as a function 
of energy and impact parameter b of the collision - these are after all often variables 
controlled experimentally. While the transformation to impact parameter representation 
suffers from some complexities, our main focus will be on collisions in the ultrahigh-enegy 
limit, E > M D , where M®~ 2 = (2ti) d - 4 /(87tG d ) gives the D-dimensional Planck mass. 
There, for many purposes, we expect the classical relation 

Z ~ £76/2 , (3.1) 

which should approximately hold more generally, to serve as a useful guide to the physics, 
though we expect precise statements to be more easily made in terms of the conserved 
quantities E and I. 

7 




Fig. 1: Scattering regimes in an impact-parameter picture; the question 
marks denote possible model dependence discussed in section 3.3. 

Fig. 1 illustrates some of the regimes that we expect to be relevant for ultrahigh-energy 
scattering, in terms of energy and impact parameter. We will particularly focus on the 
Born regime, the eikonal regime, and the strong gravity, or "black hole" regime. 

3.1. Born and eikonal 

The best-understood regime is the Born regime, corresponding to large impact param- 
eters/small angles. Here, the elastic scattering amplitude, corresponding to single graviton 
exchange, has been given in ( |2.5|) ; one may also consider corrections due to soft graviton 
emission 

As the impact parameter decreases, or the energy increases, diagrams involving ex- 
change of more gravitons become important. The leading contributions at large impact 
parameter are the ladder and crossed ladder diagrams, which can be summed to give the 
eikonal approximation to the amplitude 0,|^,[JT]|52|] .i This can be written in terms of 



2 One may inquire about UV divergences of loop diagrams. However, these are short distance 
effects, for which we assume there is some UV regulation; for example, string theory might serve 



this purpose, or even supergravity, if it is perturbatively finite [32]. 
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the eikonal phase, which arises from a Fourier transformation converting the tree-level 
amplitude into a function of a variable naturally identified as the impact parameter: 

a n (3-2) 

4tt G d s k ' 



(D - 4)^_ 3 



x 



.D-4 ' 
1 



where q± is the transverse momentum transfer and where 

27r (n+l)/2 

Qn = r[(n + l)/2] (3 - 3) 
is the volume of the unit n-sphere. The eikonal approximation to the amplitude is then 

iT eik (s,t) = 2s J d D - 2 x ± e~ ici± -^(e i ^ x± ^ - 1) , (3.4) 

expressing the amplitude in an impact-parameter form. From (|3.4| ) , one sees where eikonal 
corrections to the Born amplitude become important, namely when the eikonal phase x 
becomes of order one. Indeed, showed that at the corresponding point via (PTT|) , the 
partial wave phase shifts become of order unity, and thus the eikonal amplitudes unitarize 
the amplitudes of the Born approximation. (Contributions due to soft graviton emission 



were also estimated in [[22].) In terms of impact parameter, this transition region is given 
by 

6~ (GdE 2 ) 7 ^ 1 , (3.5) 

as is illustrated in Fig. 1. It is alternatively described as the region where the momentum 
transfer is of order the inverse impact parameter, 

V^i ~ i . (3.6) 

In general, eikonal approxmations are expected to capture semiclassical physics. In 
the high-energy gravitational context, the semiclassical geometry is the collision of two 
Aichelburg-Sexl shock waves, and various evidence supports the correspondence between 
( |3~4| ) and this picture [^,0-Q . In particular, the saddle point of ( |3.4| ) gives a classical scat- 
tering angle 



a 1 9 

E do 



mi-:) 1 ' 



(3.7) 



b 

matching that of a test particle scattering in the Aichelburg-Sexl geometry. Here, we have 
introduced the Schwarzschild radius corresponding to the CM energy, 

1 /kcE^ 1 '^ 



where „ . 

2(27r)^- 4 , . 
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Fig. 2: The H-diagram, which provides a leading correction to the eikonal 
amplitudes as scattering angles approach 6 ~ 1. 



One finds that corrections to the ladder series become important when y^i ~ E, 
or alternatively when the scattering angle reaches 9 ~ 1 . Eq. ( |3.71) shows that this happens 
at impact parameter comparable to the Schwarzschild radius, b ~ R(E), as pictured in 
Fig. 1. A schematic argument for this follows from power-counting. Consider a diagram 
arising from a graviton tree attached to the external lines. Each graviton vertex gives a 



factor yGr). Those connecting to external lines are accompanied by a ^fs. The remaining 
dimensions come from internal (loop) momenta. For the processes in question, these have 
characteristic value! k ~ 1/6. This counting then produces a power series in (R/b) D ~ 3 . A 
leading such correction, the H-diagram, which has been discussed in is illustrated in 
Fig. 2. One can alternatively understand this expansion by thinking of the external lines as 
classical sources; using standard power-counting techniques ||34|, one can easily show that 



the H-diagram is <D[(G D E) 2 /r 2 ( D -V} compared to one graviton exchange, if the distance 
between the sources is r [^5|. Using GdE ~ R D ~ 3 and taking r ~ b then yields the same 
expansion parameter. In terms of the semiclassical geometry, at impact parameters b ~ R, 
one forms a trapped surface [^,^7|, and hence a black hole. 



3.2. Strong Gravity 

Since corrections to the eikonal amplitudes give terms that differ from the eikonal 
amplitudes by powers of [R(E)/b] D ~ 3 , the region where a classical black hole forms appar- 
ently corresponds to a manifest breakdown of the perturbative expansion; it is not even 



3 Indeed, in the eikonal regime, the dominant term in the exponential series of (3^) occurs at 
order iV ~ GDs/b D ~ 4 , corresponding to a characteristic momentum k ~ y/^t/N ~ 1/6 in each 
internal line of the ./V — 1-loop Feynman ladder diagram. 
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asymptotic. We can also parameterize this in terms of a critical angular momentum, given 
by 

/ ~ L(E) = ER(E)/2 . (3.10) 

One might be tempted to believe that a quantum treatment of the evolution can still 
be given by performing an expansion in fluctuations about a shifted background - that 
of the semiclassical black hole. However, the problem of the singularity guarantees this is 
not a complete description. Moreover, even evolution on spatial "nice slices" that avoid 
the singularity is problematical, given that a standard field theory treatment of it leads 
to the information paradox.^ This suggests that the boundary of this regime represents a 
correspondence boundary, analogous to that for example between classical and quantum 
mechanics, beyond which local quantum field theory does not give a complete description 
of the dynamics ||40|| . In particular, the unitary evolution which we are assuming, in which 



the quantum information must escape the black hole while it is still comparable to its orig- 
inal radius[^TJ, suggests that the nonperturbative dynamics unitarizing the physics is not 
local with respect to the semiclassical geometry - a sort of "nonlocality principle j23| , p6| ■ " 
(This then fits with the proposed parameterization of part of the correspondence bound- 
ary given by the locality bound ||25| , |IT| , |2B|| : namely local field theory fails for multi-particle 
states whose wavefunctions are concentrated inside a radius of size R(E), where E is their 
combined CM energy.) 

While we do not have the means to calculate quantum amplitudes in this regime,!! we 
can infer some of their properties if we believe that the semiclassical picture of formation 
of a black hole and its subsequent evaporation provides a good approximate description 
of the physics when addressing certain coarse-grained questions. Specifically, ref. [p2 



parameterized certain features of the corresponding S-matrix, and we will improve on the 
corresponding "black hole ansatz" in subsequent sections. 



For reviews, see p8| , |39| ] . 

5 Ref. [13 1 has suggested analytic continuation of the perturbative sum giving the amplitude in 



the region b > R. However, one might at best expect such a sum to approximately reconstruct the 



semiclassical geometry, as in [42]. Then, in particular, it is not clear how the resulting prescription 
would give unitary amplitudes that escape the usual reasoning behind the information paradox, 
which as we have summarized, apparently requires new dynamical ingredients. Indeed, this paper 
elaborates on the view that local QFT cannot fully capture the physics of the strong gravitational 
regime semiclassically associated with black hole formation. 
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Of course, investigating the internal dynamics seen, e.g. by observers falling into black 
holes, and reconciling that with outside observations such as described by an S-matrix, re- 
mains a challenging problem. Ref. has argued for flaws in the "nice slice argument" for 
information loss, of two origins. First, attempts to measure the nice slice state at a level 
of precision appropriate to investigate information loss lead to large backreaction on the 
state. Secondly, fluctuations e.g. in the Hawking radiation are argued to lead to fluctua- 
tions in the nice-slice state after long times. We expect that sharper investigations should 
follow from use of proto-local observablespBf, but ultimately the full non-perturbative dy- 
namics of gravity is plausibly necessary in order to give both a complete picture of infalling 
observers and of reconciliation of their observations with a unitary S-matrix. 

3.3. Other regimes 

Before turning to further description of the strong gravity regime, it is important to 
note that at impact parameters larger than b ~ R(E), other features of the dynamics can 
become relevant. Indeed, some have argued that this indicates other dynamics besides 
strong gravity is a dominant feature of high-energy scattering. To give an example, in the 
context of string theory, with string mass M 8 t, it is possible to make long strings with 
length / ~ E/Mg t . In fact, such processes are highly suppressed, but pointed out that 
such amplitudes receive other important string corrections through "diffractive excitation" 
beginning at impact parameters of size bt ~ M- x {E/M st ) 2/( < D - 2 \ Indeed, g§ proposed 
that this effect may provide important corrections to a picture of black hole formation; if 
true, this would likely obscure a strong-gravity interpretation of the regime b < R(E). 

Refs. [|TT|JT^1 investigated these effects more closely. Indeed, as pointed out in [fTTfl , 



a simple picture of the origin of these effects is string excitation arising from the tidal 
impulse of the gravitational field of the other colliding string. Moreover, [[12] investigated 
the evolution of the corresponding string states. For impact parameters b t ^> b ^> R(E), 
the asymptotic state of the string is indeed highly excited as a result of this tidal string de- 
formation. However, for impact parameters b < R(E), the timescales of horizon formation 
and string excitation differ significantly. Roughly, in a semiclassical picture the trapped 
surface forms before the tidal excitation causes significant extension of the string. Thus, 
one seemingly produces a configuration described as a pair of excited strings inside a black 
hole; in this context there is no clear reason to believe that string extendedness would lead 
to significant modification of the black hole description of the dynamics. Likewise, there is 
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not a clear mechanism for string effects to provide the necessary nonlocality with respect 
to the semiclassical picture, to allow information escape. 

Indeed, one can imagine a similar dynamics being relevant for collisions of other 
composite objects - hydrogen atoms, protons, etc. Specifically, when tidal forces reach a 
size sufficient to excite the internal degrees of freedom of the object, asymptotic states will 
be excited states. Thus, there can be model-dependent tidal excitation effects. However, 
once impact parameters reach the regime b < R(E) (and for sufficiently large E), such 
effects are not expected to prevent black hole formation. Since these model-dependent 
tidal-excitation effects do not appear to contribute fundamental features to the story, we 
will largely ignore them in the following discussion. 

Another regime that has been of much interest in string theory discussions is that near 
the string energy, E ~ M s t, where one might expect to initially see weakly-coupled string 
excitations. This region lies in the lower left corner of Fig. 1. One expects such excitations 
to merge into black holes at a "correspondence point|3^]" where R(E C ) ~ 1/M st . Our 
focus will be on higher energies. 



4. The strong gravitational regime 

We currently lack a complete quantum description of the strong gravitational regime. 
However, we will assume that the quantum description of this regime must be compat- 
ible with certain features following from a semiclassical picture of black hole formation. 
If one accepts such a viewpoint, and moreover assumes that the microphysical evolu- 
tion is unitary, these combined assumptions potentially provide interesting constraints on 
the dynamics - particularly in view of the preceding statements that unitary evolution 
is apparently incompatible with evolution that is local with respect to the semiclassical 
geometry. 

4-1. Black hole formation 

We begin by recalling basic features of black hole formation in a high-energy collision, 
which has been extensively studied as a phenomenological feature of models with a low 
Planck scale US© 

Consider a high-energy collision of two particles, with CM energy E ^> Mjj. Let us 
moreover assume that the wavefunctions of these particles are gaussian wavepackets with 



For a review with some further references, see [48|. 
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characteristic size Ax, and that these collide with an impact parameter b < R(E); for 
large E, we may take Ax <C R(E). 

In the classical description of this process, a trapped surface will form in the 
geometry [ ^9j|36| , signaling formation of a black hole, and as a result of the small cur- 
vatures, one expects a corresponding statement in a semiclassical approximation to the 
quantum dynamics |$7|. Not all of the collision energy is trapped in the black hole, which 



is initially rather asymmetrical, and radiation (soft gravitons, gauge fields, etc.) will escape 
to infinity during the "balding" process in which it settles down to a Kerr black hoi J of 
mass M. The time scale for balding is of order Tf orm ~ R(E), and for impact parameters 
sufficiently below R(E), the amount of energy lost is an 0(1) fraction, but not large (e.g. 
< 40%), thus M « E. 

Subsequently, the black hole will radiate, initially preferentially radiating states that 
lower its spin. The characteristic energy of radiated particles is the Hawking temperature, 
Tu ~ 1/R{M), and roughly one quantum is emitted per time r ~ R(M). 

4-2. Black holes as resonances 

We will thus think of the black holes that form after r/ orm as resonances ||22|| . Since 
the width for such a state to decay (typically into a lower-energy black hole) is F(M) ~ 
1/R(M), this is a limit to the sharpness with which we can define the energy of the black 
hole. However, black holes with M ^> Md are sharp resonances in the sense that 

T 1 1 

— ~ ~ <C 1 (41) 

M RM S(M) ' y ' J 

where S(M) is the Bekenstein-Hawking entropy. 

We will assume that the number of possible black hole resonances is given by this 
entropy. To be more precise, let us assume that the number of black hole microstates with 
energies in a range (M, M + AM) is 

AJV(M) = B(M)e s{M) R(M)AM , (4.2) 

where B(M) is a possible prefactor that is dimensionless and is expected to have much 
more slowly-varying energy dependence than the exponential. Thus the density of black 
hole states is of the form 

p(M) = RBe s{M) , (4.3) 



In models with gauge charges not carried by light particles, the black hole can also carry 
charge. 
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and the total number of states with energy < M is J\f(M) ~ B (M) exp{S (M)} . The 
spacing between the states is clearly much smaller than their widths. Let us label the 
states in the interval (M, M + 1/R) as 

|M;/), (4.4) 

where I = 1, • • • , AJ\f(M) ~ exp{S'(M)}. We may further refine the description to project 
on angular momentum eigenstates, with angular momenta /. In that case, the entropy 
entering the preceding formulas is expected to be 

= (4.5) 



where R(E,l) is given bypD[ 



AM 2 

For small I, this gives an expansion of the form 



16nG D M 
(D-2)Q D - 2 ' 



(4.6) 



S(M, I) = S(M, 0) ^1 - const. ^\ . (4.7) 

4-3. Black hole spectrum and evolution 

Let us explore in more detail the quantum states formed in a collision, which could 
be either a two-particle collision with a CM energy « E, or an n-body collision. Note 
that one can also form a black hole of mass M by producing a higher-mass black hole in 
a collision with E ^> M, and then waiting for that black hole to evaporate to M. 

Consider general initial multi-particle (but not black hole) states; these can be labeled 
by energy, momentum, generalized partial waves, and asymptotic species and spin content. 
Let us work in the CM frame, and ignore the effects of particle spin. Some subset of the 
states, denoted \E;a)i n , will form a black hole; examples are the two-particle states de- 
scribed above, which classically do so, and thus are expected to have probability essentially 
unity for black hole formation. 

This means that a state! \E\a)i n can be rewritten in terms of states that at a time 
just after formation corresponds to a combined state of black hole and balding radiation; 
let us choose an orthonormal basis \E';i) ra d for the latter, and thus write 

\E; a) in = A ( E > M )an\M; I)\E - M; i) rad ; (4.8) 



A more careful treatment uses narrow wavepackets. 
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here we neglect the possibility of a small component on states that are not black holes. 
In principle we can project on a definite state of the radiation, yielding a pure black hole 
state: 

rad (E-M;i\E;a) in = ^ A(E, M) a n\M\ I) . (4.9) 

I 

In a generic black hole basis we expect the amplitudes A(E, M) a n to be of order 
e -S(M)/2^ corresponding to the fact that from (O) we expect there to be 0(e s ) states. 



The space of states in ([4.9|) can be combined to form an orthonormal basis for a subspace 
of black hole states, denoted \M; A), and labeled by the initial and radiation state labels. 
However, this basis will not span the space of all black hole states, since Q4.9| ) yields too 
few states. Indeed, note that there are arguments (extending |M]]) that only of order 

(g-2)(Q-l) 

exp{E D(c- 3 ) } (4.10) 

states can be formed from collapse of matter of energy E; thus a should have such a 
range. If one also accounts for the balding radiation, as above, there are more states 
that can be accessed through their entanglement with this radiation. Typical radiated 
quanta have energies ~ 1/R, and given the radiated energy E — M, this yields an entropy 
~ R(E)(E — M) oc .E^- 0-2 )/^- 0-3 ). This exponentiates to give the number of states over 
which the index i can range. However, this is still far fewer than the exp S(M) black 
hole states, since typically M > E/2. Thus, the number of states that are "accessible" in 
the collision at energy E is far less than the number of possible states of the black hole. 
We can label a basis for the remaining complementary black hole state space as \M; A). 
One expects that one approach to accessing these states is to form a black hole of mass 
M' > M in a higher energy collision, and then allow it to evaporate down to mass M. In 
doing so, internal states of the black hole become entangled with the state of the Hawking 
radiation, like in the preceding discussion of balding radiation.! For large enough M', this 
gives e s ( M ) independent accessible states. For many purposes, it is simplest to forget the 
balding radiation, which as we have explained does not appear to play a particular central 
role, and in a slight abuse of notation, think of the labels A as corresponding to the initial 
states from which the black hole formed. 

We can likewise label the possible n-body out states, representing the final decay 
products of a given black hole, as \E, a) out . In a similar spirit to the preceding discussion, 



9 One can in principle "purify" such states by projection on definite states of the Hawking 
radiation, as with the preceding projection of balding radiation. 
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we could choose a basis of black hole states labelled by this out-state description. Again, 
we expect the matrix elements between the preceding basis and this one to generically 
have size exp{— S(M)/2}. Correspondingly, the amplitude for a given initial black hole 
state to decay into a given final state of the Hawking radiation will be of generic size 



out ( 



(M,a\M,I)\ ~e~ s(M)/2 . (4.11) 



The quantum description of black holes as a decaying multi-state system has analogies 
to other such systems, like KO — KO mesons. In the assumed unitary dynamics, an initial 
black hole state \M; I) can both mix with other states with the same energy, and with 
states that are in the continuum, which consist of a lighter black hole together with radiated 



quanta. One might expect, via a Weisskopf-Wigner ||52|| approximation, that evolution 
in the Hilbert space of black hole states with mass ~ M is governed by an effective 
Hamiltonian: 

Am;I) = H\M;I). (4.12) 
at 

Though conceivably more general dynamics is needed,0 this exhibits possible features of 
black hole evolution. Due to the decay, the hamiltonian is not hermitian in this subspace, 
and in general takes the form 

H I j = M I j- % -T IJ , (4.13) 
where Mu and Yjj are hermitian matrices. In general, these will not commute. 

4-4- Exclusive processes 

If one considers in particular an exclusive process with two-particle initial and final 
states \pi,P2)in, \P3,Pi)out, such as pictured in Fig. 3, one thus expects the intermediate 
black hole states to contribute to the S-matrix as 

out(ps,P4\pi,P2hn = (2n) D S D (J2p) X>3,P4|J> (^) (I\Pl,P 2 ). (4-14) 

IJ \ ' JI 

(Note that in the bases adapted to in or out states, described in the preceding subsection, 
the indices are expected to only range over ~ S(E) values.) If Mjj and Tjj do not 
commute, Hjj cannot be diagonalized by a unitary transformation, but we will assume it 



10 In particular, we don't expect H to necessarily be a hamiltonian constructed from a local 
lagrangian. 
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Fig. 3: Schematic of a black hole as a resonance in 2 — > 2 scattering. 

can be diagonalized by a more general linear transformation. The eigenstates \M;I) are 
then not orthogonal; 

(M;I\M;J)=g n (4.15) 
for some gjj ^ djj. In such a basis ( |4. 14f) becomes0 



out 



E-Hf 
ij 



where Hj = Mj — iTj/2 are eigenvalues. This will produce a sum of terms of Breit- Wigner 
form contributing to the amplitude. However, the sum itself will not, in general, take the 
Breit- Wigner form. 

In the case where the particles being collided are the narrowly- fo cussed wavepackets 
that we have described, one plausibly expects the corresponding amplitude to be of size 



out ( 



<a|6) in |~e-^)/ 2 . (4.17) 



The reason for this is that for such wavepackets the amplitude to form a black hole is 
essentially unity, and the amplitude for it to decay back to a two-particle state is of size 
given by ( 4.11 ). Note that our discussion suggests a resolution to questions raised [[KJ about 



the relation of intermediate black holes to Breit- Wigner behavior. One has 0(1) amplitude 



11 The form of this equation may alternately be simplified through the definition of a dual basis, 
(I d \=9j}{J\. 
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to form some black hole state; in a generic basis for black hole states, this is a superposition 
with 0(e~ s / 2 ) coefficients, although, as indicated in the preceding subsection, one can 
choose a special basis where black hole states are labeled by the initial states that created 
them. Thus, the amplitude to form a generic black hole state from a two-particle state is 
~ e -5 / 2 , as is the amplitude for a generic black hole state to decay back into a two-particle 
state. 

One might ask whether there could be any larger contributions to the 2 — > 2 am- 
plitude, due to processes that avoid black hole formation. For example, our gaussian 
wavepackets will have tails at large impact parameter. However, these have probability 
of size exp{ — {R/ Ax) 2 } at b ~ R. The width Ax is constrained by Ax > 1/E, but this 
constraint produces a quantity merely of size > exp{— S 2 }. 

While we can't at present rule out other such effects, none have been identified. An- 
other test of this statement comes from scattering of a particle of high energy E off a 
preexisting black hole in the relevant range b <C R; here the amplitude 1Z for reflection is 
also exponentially suppressed [(53fl: 



U ~ e~^ ER . (4.18) 



It is thus plausible that the amplitude for the classically predzcted [p9| , |3q , p7| black hole 
formation process only receives corrections that are exponentially suppressed at least to 
the level ( pTj) . 



5. Partial wave amplitudes 

In this section we restrict attention to 2 — > 2 scattering, in a partial wave basis, and 
investigate consequences of the preceding picture and related considerations. For simplicity, 
we focus on scattering of one species of spinless particles. The initial two-particle states 
will be labeled by just their energy and angular momentum /, and the scattering amplitude 
is of form 

S t (E) = e -Wi(E)+2i5UE) (51) 
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5.1. Strong gravitational regime 

As outlined above, for impact parameters b <C R(E), or correspondingly angular 
momenta Z <C L(E), the amplitude for such a state to form a black hole with total angular 
momentum Zbh ~ I is expected to be of order unity. 
Absorption 

In the 2 — > 2 process that goes through the black hole channel, Zbh = I- From ( |4.11|) , 
we note that the amplitude for the given resonance \E, Z) to decay back to a two-particle 
state is ~ e -s(E,i)/2_ 

As in the preceding section, it is plausible that processes avoiding black hole formation 
in the regime Z <C L are exponentially suppressed at least to this level. Arguments for that 
build on the preceding ones, together with the properties of partial-wave wavepackets. 

For example, consider a wavepacket with definite angular momentum in the relative 
coordinates between the two particles: 

1>im{x) = J dE ^±^- e- iEt Y lm (n)f(E), (5.2) 

where Ji+\ is a Bessel function, Y/ m (f2) are D — 2 dimensional spherical harmonics, and 
f(E) is a gaussian wavepacket with width AE. Asymptotics of Bessel functions for large 
order and argument (see eq. 8.41.4 of ||54| ) then show that for Z <C Er, 



J,-x(Er)^ \/^cos 



2 4 



(5.3) 



with subleading corrections consisting of terms suppressed by powers of (Z + A) / Er times 
sine or cosine functions of the same form. Thus combining ( |5.2| ) and ( |5.3| ) gives a 
wavepacket that is gaussian in t ± r with width Ar ~ 1/ AE, and subleading terms are 
similarly gaussian. 

A related argument comes from the relation between the partial wave representation 
and impact parameter representation |55[]. Specifically, if f(b, s) is the amplitude in impact 



parameter representation, then at high-energies one finds the corresponding partial wave 
amplitude [j56l , |5"7H 



m = m , + X)/E,s) + ^^l {2v+xWE + ..., (5.4) 
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where A is a numerical coefficient, indicating that in the high-energy limit, localization in 
angular momentum corresponds to localization in impact parameter, as expected.0 

A final argument comes from the behavior of partial waves scattering from a preex- 



isting black hole; |53[ argues that their reflection amplitude in the limit ER 3> 1 is of size 

(EH- 

Based on these, and on the discussion of section four, we thus conjecture that in the 
regime / <C L(E), the 2 — > 2 amplitude is indeed exponentially small in the entropy, and 
arises mainly due to such a strong gravity channel. These statements suggest additional 
rationale for the "black hole ansatz" of ||22|| , that in this regime 

\Si(E)\ =e" 2A ~exp{S{E,l)/2}. (5.5) 

Notice that this behavior has two characteristic features. The first is the exponential 
strength of the absorption. The second is the long range of the absorption, which is 
characterized by the growth of L(E) with energy. Even should the preceding arguments 
regarding the strength of the exponential suppression be evaded, we expect the feature of 
significant absorption at long range to persist. 
Phase shifts 

We have suggested that the amplitude is essentially unity for a given initial two- 
particle state with / <C L(E) to enter the strong gravitational regime. In 2 — > 2 scattering, 
one might therefore expect that in each energy range (E,E + 1/R) we form one of the 
AJ\f(E } I) black hole states^! with the corresponding energy and angular momentum. This 
would correspond to a density of "accessible" states 

Pacc {E,l)nR(E) . (5.6) 

(This value would be less relevant for 2 — > N scattering, where, as we have argued, more 
states may be accessible and entangle with the balding radiation.) Notice that this would 
imply that the total number of such accessible black hole states of angular momentum I 
and energy < E is given by 

Kcc(E) = [ pacc(£, l)dE « S(E, I) . (5.7) 
Jo 



12 The series ( |5.4D may be regulated by considering incoming wavepackets instead of plane 
s. 

As noted, this state is a superposition of states of a generic basis with coefficients of size 



waves. 
13 



0{e 



S/2s 



21 



Consider the parametrization (|4.14| ) of the contributions of intermediate black hole 
states. If the matrix Hjj were diagonal in the "in"-state basis \M; A), discussed in section 
four, then we would expect a contribution to the amplitude of Breit-Wigner form: 

where Sj, is a "background" value. Then, the phase 5i(E) would increase by tc as we pass 
through each such "accessible" (or strongly coupled) resonance, and correspondingly, the 
combined effect of resonances at increasing energies would give 

Sf ag (E) = ttA/; cc (£, « ttS(E, /), (5.9) 

as with Levinson's theorem for single-channel scattering. Note also that such a result 
would yield a decay time ddi/dE ~ R(E), compatible with the width T ~ 1/R. 

However, we see no reason to expect Hjj to be diagonal, and so consider phase shifts 
of a more general form, which we parameterize as 



6i(E) =irk(E,l)S(E,l) (5.10) 

where k(E, I) varies more weakly with energy than S(E, I). One might expect k(E, I) > 
(corresponding to time delay) due to the attractive nature of gravity. Indeed, in scatter- 
ing off a pre-existing black hole the gravitational field introduces a positive phase shift 
relative to scattering from the angular momentum barrier. We will investigate additional 
constraints on k(EJ) in subsequent sections. 

To summarize, combining (|5.9|), ( |5.10|) suggests that the partial wave amplitudes in 
the strong gravity regime take the form 

/f G ( S )~|{l-e X p^-^(E,0[l-4 7 rzA ; (E,0])} • (5.11) 

Notice that this expression differs from that of [22] ; that analysis did not take into account 
the role of inelasticity and accessibility of resonance channels. Thus ( |5.11| ) comprises an 
improvement of the black hole ansatz of [p2~ |. 
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5.2. Born and eikonal 



One can likewise infer properties of the partial waves in the longer-distance regimes, 
where the Born or eikonal approximations are expected to be valid. In particular, ref. |]22| 
computed the eikonal phase shift, 

ik ^(D-2)T[(D-4)/2]L(E) D - 3 E D ~* 
61 {E) = 8T[(D - l)/2] T®- 1 I® 11 * ' (5 ' 12) 

and checked that the eikonal amplitude unitarizes the Born amplitude, which is the leading 
term in an expansion in Si, as expected. Thus the transition from Born to eikonal regimes 
occurs in the small angle regime / ~ E^ D ~ 2 ^^ D ~ A \ Notice that the phase shifts are indeed 
positive definite, as expected from the attractive nature of gravity.0 The correspondence 
between the eikonal amplitudes and the semiclassical picture ]8|,|]-[| suggest the utility of 
the eikonal description until I ~ L. 

For decreasing impact parameter/increasing scattering angle, different effects can con- 
tribute to absorption. A generic effect is soft-graviton bremmstrahlung. This was esti- 
mated in to give a contribution of size 



L(Ef D -yi™-™ ~ fL-^ . (5.13) 



Note that this matches onto the energy dependence of (O) at I ~ L, which also fits with a 
picture where a non-negligible fraction of the collision energy can be emitted in the balding 
radiation. 

As noted in section three, there may be other less-generic effects, e.g. due to excitation 
of internal degrees of freedom of the colliding bodies. In string theory, such an effect is the 
"diffractive excitation" or "tidal string excitation" explored in P-|5|, pT| , p!2| . But, as noted, 
we do not expect such effects to prevent amplitudes from matching onto those of the strong 
gravitational regime. 

5.3. Combined pictures 



14 This is the case provided D > 4. The four dimensional case suffers from Coulomb-like 
singularities, requiring the usual inclusive amplitudes, avoided in this paper by working in higher 
dimensions. 
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Bremsstrahlung 



Black Hole Regime 



D-3 

Fig. 4: Absorption coefficients at a fixed angular momentum as a function of 
the CM energy. 

A(E) 
L,. - 



Black Hole Regime 



Bremsstrahlung 



L, 

Fig. 5: Absorption coefficients at a fixed CM energy as a function of angular 
momentum, with L c = L(E). 

We can thus suggest combined pictures describing the weak and strong coupling 
regimes. The results (|5.13|) and ( |OD suggest energy and angular momentum dependences 
of the absorptive coefficients (3i as pictured in Fig. 4, Fig. 5. 

While the phase shift is well-studied in the eikonal regime, as we have indicated, 
we have less information in the strong gravity regime, but expect an increase bounded by 

5i(E) ~ £(-D-2)/(.D-3) 

as in (|5.10|) . Sketches of energy and angular momentum dependence 
are given in Fig. 6, Fig. 7. 
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Fig. 6: Phase shift for fixed angular momentum as a function of the CM 
energy. 
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Fig. 7: Phase shift for a fixed CM energy as a function of angular momentum, 
with L c = L(E). 

6. Momentum space amplitudes 

We now ask what properties of momentum space amplitudes can be inferred from the 
preceding discussion. In section two, we noted the collapse of the Lehmann ellipse, and 
in particular that convergence of the partial wave expansion cannot extend past t = 
to positive t. Likewise, continuation of s to complex values with fixed real t < would 
correspond to complex cos 9, outside the convergence region. These and related limitations 
restrict our ability to prove results that follow in massive theories. However, we have argued 
that the expression for the partial wave coefficients, ( |2.10| ), is expected to be well-defined 
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and finite. This means that properties of the fi(s) are those of the corresponding integral, 
and this in turn constrains the behavior of T(s,t). 

Additional information about the momentum space amplitudes comes directly from 
their eikonal approximation, ( |3.4j ). At very small angles, this expression reduces to the 
Born amplitude, ( |2.5p . The match between the Born and eikonal regimes occurs near 
X ~ lj corresponding to t ~ — s _2 /(- D_4 ) or 

&B/E ~ £ (D _ 2)/(D _ 4) - 

The asymptotics of the eikonal amplitude at larger angles follows from performing the 
integral over angles in ( p.4j) , which yields 

POO 

iT eik (s,t) = -2z S (2tt)^- 2 )/ 2 ^ (D - 4)/2 / dxxxf - 2)/a J i n =& {q±x ± ){e i ^^ - 1). 

Jo ~ 2 

(6.2) 

Then, combining the Bessel function asymptotics (5.3) with a saddle-point approximation 
of the integral gives an asymptotic amplitude of the form 

T eik ~exp{^[ S (-t)^- 4 )/ 2 ] 1/(D - 3) } • (6.3) 

This exhibits some interesting features - such as nonpolynomiality - that we will return 
to in the next section. 

One may also inquire about implications for T of the strong gravity behavior outlined 
in the preceding section. Recall that the physical features of that behavior were 1) sig- 
nificant scattering, and moreover absorption, to an angular momentum that grows with 
energy as I ~ L(E), 2) strong absorption for large E and I <C L(E), and 3) potentially 
rapid growth in the phase, (|5.10|) . 

For 5i = Pi = 0, ( |2.13| ) gives // = 0, so the first feature implies nonvanishing fx to 
I ~ L(E); significant absorption moreover implies that fx ~ i/2. These become conditions 
on the integral 

(« D «,^ T MM)] = |f, (6.4) 

where t(s, 9) is given by ( |2.9| ). However, a direct statement about T in the strong gravity 
regime s ~ — t, is not easily inferred from the significance of the right side of ( |6.4j ), since 
the integral in particular receives a contribution from the Born regime. For 9 < O^/e an d 
/ < L, one has Wcl and can use the small-angle approximation 

Ct { l-eV2),Ctm(l-10^) . (6.5) 
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The Born contribution to (6.41) is thus of size 



} B/E T?2 1 

^ D_ V ~ f*=*- (6 - 6) 

This shows that one expects a contribution to partial wave amplitudes from both the Born 
and eikonal regions that is significant at angular momenta I < L(E). 

Indeed, a related fact is that the cross- section due to this small-angle scattering is 
expected to be large as compared to that of the strong gravity region, 

a SG « tt[R(E)} d - 2 ~ eP-W-Q. (6.7) 

For D > 6, where the small angle contribution converges, it can be estimated using the 
impact parameter where Born and eikonal match, giving @ 

a B/E ^E^ D - 2 ^^. (6.8) 

Large growth of j3i and 5i with energy imply that fi — i/2, or dfi/ds, are small, 
and rapidly oscillating. Eq. (|6.4| ) thus indicates that T(s,t) correspondingly has rapid 
falloff and oscillations. Moreover, we see that exponential falloff of fi — i/2 would indicate 
precise cancellations between the contributions of T(s,t(s,6)) in the Born, eikonal, and 
strong gravity regimes; as we have discussed, physical aspects of the scattering such as the 
analogy with scattering from a fixed black hole suggest such falloff. 

A sharper statement arises if one considers continuation of (|5.11| ) into the complex s 
plane. This form for fi(s) suggests that generically it would grow exponentially somewhere 
in the complex s (or E) plane. In particular, for small enough k, one finds exponential 
growth in the s upper half plane (UHP) < Args < it: for constant k, this would occur 
for 

k < — tan — — — -, (6.9) 
4tt 2(D-3) K ! 



and likewise for the example of a decreasing power, k oc E~ p . By fl6.4| ), this corresponds 
to exponential, thus not polymomially bounded, growth in T(s, t) for complex t. While 
with the specific functional form ( |5.11|) , a phase that is too small leads to growth that is 
not polynomially bounded, it is conceivable that a more complicated analytic structure of 
the exact amplitude avoids this conclusionJil 



15 Though, with added assumptions like hermitian analyticity /dispersion relations, one may 
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7. Analyticity and crossing 

We have investigated aspects of unitarity, particularly via the partial wave expansion; 
we now turn to analyticity and crossing. 

Consider scattering of two massive particles of mass m coupled to gravity. We might 
imagine these to be an e + e~ pair, although to avoid complications of spin we will treat 
the scalar case. Another specific context to contemplate, if in a string theory context, is 
scattering of a DO — DO pair. 

First, consider behavior for fixed real t < 0, as a function of s. The two-particle cut 
in the s-channel begins at s = 4m 2 . However, one can also have such a pair annihilate to 
two or more gravitons (in the absence of a net conserved charge), implying multiple cuts 
beginning at s = 0.0 Likewise, there are multiple u-channel cuts beginning at u = 0. 
Given 

s + t + u = Am 2 , (7.1) 
we find that the u-channel cuts, for fixed t, originate at 

s = 4m 2 -t, (7.2) 

and are taken to extend along the negative s axis. Thus, these cuts overlap those from 
s = - there are branch cuts running all along the real s axis, with no gap between them, 
unlike the massive case. These features of massless theories weaken some of the constraints 
present in massive theories. 

We likewise expect singular behavior at t = 0; we have noted the Coulomb pole 
there, but one might find a more general singularity (e.g. branch point) when higher- 
order processes are accounted for. As we have already described, this prevents the usual 
continuation along the real axis from t < to t > 0, that is a useful tool in massive 
theories. 



possibly generalize methods of [58,59] to show that the exponential falloff in ( 5.1 1| ) implies a 



lower bound on the phase, e.g. 8 > logs, given a polynomial bound in the UHP; also, certain 
analyticity assumptions together with this falloff might possibly be used to prove violation of 



polynomial bounds in some region, with methods like in [ 30,61] . We leave these for future investi- 
gation. (Notice that in QFT we do not expect such a strong absorptive behavior, thus polynomial 
boundedness is expected to lead to a phase bounded above by logs.) 
16 One might also contemplate the possibility of worse behavior, e.g. ~ e _1 ^ sP for some p. 
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7.1. Crossing symmetry 



For real so > 4m 2 , the physical amplitude with s = sq, t < is assumed to arise 
from the analytic function T(s,t) with s = sq + ie in the limit e —>■ + . By the maximal 
analyticity hypothesis, T only has singularities dictated by unitarity, so can be continued 
throughout the s UHP; likewise for fixed s, one can continue in t, avoiding singularities. 

In a massive theory, at small t < 0, one can continue in s across the real axis, through 
the gap between the cuts. This allows one to define the amplitude for s = si — ie, for large 
negative real s\, which by (|7.1|) corresponds to u-channel kinematics. Crossing symmetry 
is the assumption that a single function T(s,t,u), with variables satisfying (|7.1|), defines 
amplitudes in all channels through such continuation. 

Clearly this specific continuation fails in the massless case, given the lack of a gap 
between the cuts. However, it appears possible to still obtain crossing, through use of a 
different path. 
The BEG path 



Such a path was given by Bros, Epstein, and Glaser in ||62|1 , as follows. First, begin at 
large so > 0, and hold u = uq < fixed. One can continue through the upper s-plane to 
e in SQ. Here, t will approach the positive real axis with a — ie; we can denote this as the t~ 
channel. Next, beginning at this point, keep s < fixed and continue t —>■ e~ l7r t. This is 
analogous to the preceding continuation, and takes t~ to u + - here the positive real u axis 
is approached from above. The combined path thus continues from the physical s-channel 
s + to the physical u-channel u + , permitting crossing.0 

7.2. Crossing and polynomial boundedness 



17 Note that one must also include a small path segment from (s,t,u) = (—so + ie, 4m 2 — 
uo + sq — ie, uq) to (—so, 4m 2 — uq + sq — ie, uq + ie). We assume this is permitted by sufficient 



holomorphy in this neighborhood, as in [33|, though more systematic investigation is conceivably 
warranted. 
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Analyticity and crossing constrain possible non-polynomial behavior, as we will now 
discuss; the reader may wish to refer to figure Fig. 8. This observation follows from the 
Phragmen-Lindelof Theorem: If an analytic function is bounded along two straight lines 
sustaining an angle —, e.g. \T(\s\)\ < M on the lines, and ifT(s) grows at most like e' s ' /3 
with (3 < a in any other direction, then in fact T(s) is bounded by M in the whole sector 
sustained by the two lines. 

Choose, for example, a = 1. Let us assume that the amplitude is quite weakly 
bounded, \T(s,t < 0)| < e' s l. Note that this bound is easily satisfied both by the eikonal 
behavior (|6.3| ) , and by behavior that could arise from growth of the strong gravity region, 
either from the large absorption coefficients |A( S )| ~ |s|( D_2 )/( 2 (' D_3 )) <c |s|, or the large 
range R(E) ~ E X ^ D ~^ which suggests behavior ||22|| (see the next section), T(s,t < 0) ~ 
e R(E)Vt^ Therefore, by the theorem, if we had a non-polynomial growth in the UHP, that 
would also require a non-polynomial growth in a straight line ie above the real axis from 
— oo to +0O. 

The region [0, +oo) corresponds to the s-channel amplitude. However, properties of 
the Gegenbauer polynomials combined with the optical theorem (see appendix) show that 
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ImT(s,t < 0) < ImT(s,0) ~ sctt(s) < s . (The polynomial bound at t = is directly 
connected to existence of a forward dispersion relation |22|], following from causality, to be 



discussed in the next section.) Moreover, we have the high-energy expression 

J d6 sim 0-3 #|ReT| 2 oc J dfl D - 2 \ReT\ 2 < J dfi D _ 2 |T| 2 oc s 3 - D/2 a 2 ^ 2 < s 3 - D/2 a T , 

(7.3) 

where proportionality is modulo numerical coefficients, and therefore the real part of the 
amplitude also must be polynomially bounded, provided it is sufficiently smooth. (Recall 
that in the strong gravitational regime the real part of the amplitude is indeed subdominant 
due to strong absorption). 

In massive theories, the (— oo,0] region is related to the u channel amplitude by 
complex conjugation.@ This follows from the property of hermitian analyticity or extended 
unitarity, which is the requirement T(s*,t*) = T(s,t)*. Notice that this implies fi(s*) = 
fi(s)* for the partial wave coefficients. If we work at negative values of transfer momentum, 
e.g. t < 0, hermitian analyticity also connects the discontinuity across the cuts due to 
threshold singularities to the imaginary part of the amplitude by 

DiscT(s, t) = 2ilmT(s + ie, t) . (7.4) 

With a mass gap, hermitian analyticity follows from reality of the amplitude below thresh- 
old, along with the Schwarz reflection property. In massless theories the status of hermitian 
analyicity remains unclear, although it seems to hold at any order in perturbation theory. 
If hermitian analyticity holds in gravity, it thus also forbids non-polynomial growth along 
(— oo, 0], and so by the above theorem, in the UHP of s. 

A conservative conjecture is that gravity respects both crossing symmetry and hermi- 
tian analyticity, and that amplitudes thus satisfy such a polynomial bound. We can check 
this in the asymptotics of the eikonal, (|673|), which does so for D > 4, as does the preceding 
strong gravity expression. 

Nonpolynomiality of amplitudes is however generally expected to give unbounded be- 
havior in other regions of s,t, and u. Indeed, one can directly see indications for such 
behavior given the partial wave coefficients (|5.11| ). For example, if k(E,l) ~ E~ p for 
some p > 0, then the strong-gravity //'s given by (|5.11|) will have polynomially- unbounded 



18 A rough argument for this follows from the relation between the continuations s — > — s and 
E — > — E; the latter corresponds to taking the complex conjugate of the amplitude. 
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behavior somewhere in the UHP Im(s) > 0. Then, ( |6.4| ) implies that T[s, t(s, 9)] must 
likewise be unbounded. Notice, though, that this is for fixed 9 rather than t; thus un- 
boundedness at large |s|e l< ^ corresponds to t ~ —\s\e l( ^. As discussed, even k(E,l) = 0(1) 
does not necessarily eliminate this behavior, though positive k - corresponding to time 
delay - decreases the region of non-bounded behavior in the UHP. Likewise, k < 0, cor- 
responding to a time advance, increases the domain of this behavior. One also observes 
unbounded behavior from the eikonal phases, ( |5.12j ). 

It is interesting that a polynomial bound in the physical region Im(s) > 0, t < (and 
correspondingly in other channels) follows from the very general assumptions that we have 
described, together with the assumption of causality in the form of the forward polynomial 
bound. We next turn to investigation of connections between polynomiality and locality. 



8. Locality vs. nonpolynomiality 

The status of locality in gravity is a very important question, given that it is one of 
the cornerstones of a local quantum field theory description of nature. Locality is also one 
of the assumptions leading to the information paradox, and conversely, certain violations 
of locality inherent to nonperturbative gravity have been proposed as the mechanism for 
information to escape an evaporating black hole [ p5|p!l| , p3| , ^0| .B 

If one is restricted to an S-matrix description of dynamics, one can ask how specifi- 
cally locality is encoded in that description. In particular, nonpolynomial behavior in the 
momenta, such as we have described, is suggestive of non-local behavior;@ a first heuristic 
for this is the observation that nonpolynomial interactions take the form e 9 in position 
space, which is clearly not local. 

For massive theories, sharper statements can be made. In particular, commutativity 
of observables outside the lightcone can be used to show that the forward amplitude is 
polynomial bounded [pTfl, |T(s,0)| < s . With a mass, such statements can be extended [|68[ 
both to t < and to complex values of t, including t > 0. 

Diffeomorphism invariance forbids local observables in gravity. It has been proposed 
that local observables are approximately recovered from certain relational protolocal ob- 
servables; initial exploration of them in effective field theory is described in ||B||J3,[rD 



19 For earlier proposals of a role for nonlocal effects, see [64,|al.35|. 

20 Although, formulations of local field theory with mild nonpolynomial behavior have been 
proposed Ml 
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However, as yet no sharp criterion for locality can be formulated in terms of these observ- 
ables, and indeed it has been argued [p5| , |43[1 that there are fundamental obstacles to such 
precise locality .0 

Nonetheless, bounds on amplitudes can also be understood from a physical perspec- 
tive, in connection with causality. This becomes particularly clear with forward scattering. 

Consider first + 1 dimensional scattering, with initial and final amplitudes related 
by an S-matrix, 

/■oo 

iP f (t) = / dt'S(t-t')^(t'). (8.1) 



Causality states that if the source vanishes for t' < 0, the response ipf does as well. In 
the complex energy plane, this arises as a result of S(E) having the appropriate analytic 
structure, and in particular the needed contour deformation arguments require that S(E) 
be polynomially bounded in the UHP for E. For example, S(E) = e~ lEr would produce 
an acausal time advance by r. 

The arguments for higher-dimensional forward scattering can be formulated in anal- 
ogous fashion; a wavepacket that scatters at zero angle should not reach infinity more 
rapidly than one that does not scatter, implying a polynomial bound, and corresponding 
dispersion relations.!! Whereas in the massive case such a bound also implies bounds for 
t 7^ 0, the collapse of the Lehmann ellipse that we have noted in the massless case obstructs 
such arguments. 

Consider, however, a physical picture of non- forward scattering, as described in e.g. 



74] ; see Fig. 9. If the scattering has a range R, a wavepacket can shorten its path by an 
amount up to R\q\/E with respect to a path going through the origin, with a corresponding 
time advance. Thus, we would expect asymptotic behavior 

S ~ e~ i ^ iR (8.2) 

which is not bounded. Note, however, that such a picture is appropriate to a repulsive 
potential. If one instead considers scattering in gravity, e.g. in the background of a high- 
energy particle, whose gravitational field is approximately Aichelburg-Sexl (see Fig. 10), 
the scattering angle is negative, and the particle receives a time delay, corresponding to 



21 For further discussion, see [71]. 

22 The relations between causality, analyticity and a well denned UV completion are interesting 
and subtle. Indeed, other strong restrictions on which IR behavior can be consistently completed 



into a causal UV theory, given existence of forward dispersion relations, are described in [72,73]. 
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Fig. 9: Illustration of scattering by a repulsive interaction of range R; the 
scattered wave at angle 6 has a path that is shorter by 2R sin | relative to 



a wave traveling unscattered through the origin, thus has a relative time ad- 
vance. 

positive phase shift, appropriate to an attractive force. If of finite range R, this corresponds 
to behavior 



In this way, long range behavior of this kind, which in the absence of a better definition 
we will also call nonlocal, does not obviously conflict with causality. The danger of a conflict 
appears even less in an attractive case which produces only time delays; correspondingly 
one has a polynomial bound for R oc E p in this case when E undergoes a small enough 
positive phase rotation. Thus, plausibly, nonlocality with time delays is consistent with the 
existence of a polynomial bound in the physical region, t < 0, Im(s) > 0. The preceding 
section also argued that crossing, hermitian analyticity, and causality imply such a bound. 




(8.3) 
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Aichelburg — Sexl Shock Wave 




Fig. 10: Illustration of scattering of a particle by the gravitational field of 
an ultrarelativistic source; the scattering angle is negative, corresponding to 
attraction, and this results in a path for the scattered wave that is longer by 
Rsin6 ~ 2R\/tu/s as compared to a wave that passes through the scattering 
center. 

While the large phase shifts and strong absorption up to large impact parameters that 
we have inferred on physical grounds might have violated such a polynomial bound in the 
physical region, we have found no evidence for such behavior. It remains possible that an 
exponential growth may emerge at fixed (real) scattering angle, other than 9 = 0. This 
however does not seem to contradict any fundamental property we know, but is another 
possible signal of nonlocal behavior.Hl 

In saying this, we should address arguments of suggesting behavior combining 
( |S.2| ) with (|8.3| ), where R = R(E), which would be naturally interpreted in terms of a time 
advance. However, this arose from a sharp cutoff in the partial wave sum and does not 



23 As noted, one might also consider the possibility, which we haven't been able to rule out, 
that amplitudes, while nonpolynomial, may have sufficiently complicated analytic structure to 
stay polynomially bounded in other regions as well. 
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account for the phase shifts. If one avoids 9 = 0, where causality requires cancellations of 
non-polynomial behavior ||22|| , we can write 



T(s,t) oc ^2(1 + A)Q A (cos6')e 2l5i(s) - 2A(s) (8.4) 

(the sum of i/2 generates a S(cos9 — 1)). Plausibly, the exact phase shifts and absorptive 
coefficients yield only time-delayed behavior, and a bound in the s UHP. 

In the preceding section, we argued that the effective range of the interaction grows 
with E; R ~ E p , with p = 1/(D — 3) for the strong gravity region, and the rough estimate 
p = 2/(D — 4), from ( |6.8|) , for the eikonal amplitudes. It is interesting to compare this 
behavior to what is commonly regarded as another indicator of unitary local behavior, the 
Froissart bound, which states 

R< R f = alogE (8.5) 

for constant a. In a massive theory, there is a direct connection between this bound and 
polynomial boundedness. Heuristically, this is seen via 

which is polynomial behavior. More sharply, the polynomial bound is used directly in the 
proof of the Froissart bound |2^,^5j. However, this proof proceeds via the partial wave 



expansion in the region t > 0, which we have argued is divergent for gravity. 

It is tempting to conjecture that there is such a direct connection between power-law 
growth of the cross section in gravity and nonpolynomiality, perhaps through appropriate 



regulation of the partial wave expansion. Indeed, as discussed in |22| and above, the 
appearance of strong absorption to L ~ i^- 0-2 )/^ -3 ) ^> ElnE implies nonpolynomial 
behavior of a truncated partial wave sum.il However, as we have argued, we expect the 
full sum to be polynomial bounded in the s UHP, even if it is not polynomial. One issue 
arising from massless modes is that we cannot neglect the tail of the partial wave expansion, 
as one does for example in theories with a mass gap, where fi decays exponentially for 



24 Note that such strong absorption directly corresponds to a cross-section with growth 
This follows from taking /J| > 1 for I <C L in ( |2.6| ) evaluated at 9 = 0; this, together with the 
large-/ asymptotics C^l) ~ l 2X - 1 /F(2X) gives T{9 = 0) ~ is {4 ~ D)/2 L D - 2 , and thus, by the 
optical theorem, fl6.7p . Of course, as we have noted, an even larger contribution to ot comes from 
the eikonal region. 
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I 3> ElogE. In our gravitational context, these large impact parameter contributions 
are central in producing the IR singularities at t = 0. Indeed, masslessness also plays an 
important role in the form of the amplitudes in the eikonal regime (where / ^> L), which 
appears to dominate the cross-section at large energies. Since the partial wave expansion 
does not converge at t > 0, the Froissart bound can be violated without collateral damage. 
We may associate this with a sort of IR/UV mixing, in the sense that the singularities in 
the IR (correspondingly the long-range character of gravity) permit a much faster growth in 
the cross section deep in the UV without conflicting with any other fundamental property. 
Notice that the eikonal amplitudes already provide us with such an example, without 
explicit reference to the strong gravity region. 

One thus finds that masslessness, and in particular singular behavior at t = 0, non- 
polynomiality, and polynomial growth of cross sections are intricately entwined. One might 
question whether all novel features follow from masslessness alone. However, given that 
one does not find power law growth R ~ E p in gauge theory, gravity appears distinctive, 
due in part to the power-law growth of its coupling with energy. One might conjecture that 
a massless theory like QED is on the borderline of locality, but gravity is in a real sense 
not local, as for example evidenced by its growth of range. Such a conjecture is certainly 
permitted without a sharper characterization of locality. 

It is interesting to consider one known approach to regulating IR behavior in grav- 
ity, namely working in an AdS background. With AdS curvature 1Z ~ /Lt 2 , the graviton 
effectively has a mass ~ \x. Correspondingly, growth of black hole radius with energy 
stops being power law once R ~ and one in particular finds evidence for Froissart- 



like behavior, R oc log-E, for scattering above this energy pql . One might likewise expect 
restoration of polynomial scattering amplitudes. However, the matter of extracting the 
S- matrix in AdS remains an open question [^0f, despite some recent progress []T8| , p!9 . 



It is very interesting that no fundamental inconsistency has yet arisen between the con- 
ditions of unitarity, analyticity, crossing symmetry, causality, and nonlocality in the sense 
described, despite the existence of nontrivial constraints arising from their combination; 
it is also moreover interesting that gravitational amplitudes could well run the gauntlet 
among these conditions. This would also been in harmony with arguments that local field 
theory breaks down in contexts described by the locality bound [p3| , |25| , p6[ , and with more 
general statements that the nonperturbative physics that unitarizes gravity (and specif- 
ically leads to unitary black hole decay) is not intrinsically local[^3|, yet retains certain 
analytic features and aspects of causality - particularly those necessary for consistency! In 
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any case, further exploration of properties of consistent quantum-mechanical amplitudes 
for gravity is certainly of great interest. 
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Appendix A. Optical theorem in D dimensions 

From the unitarity of the S'-matrix we have 

T aP -T* pa = iJ2 I (2n) D d$ N T aN T* pN (A.l) 

N J 

where we take a, (3 to be the initial and final two-body states with p a = p\ + p2^Pp = 
P3 +P4: and the sum runs over all possible A^-particle states allowed by the symmetries and 
conservation of energy and momentum. Here we use the Lorentz invariant normalization 
of states, 

(k\k') = (27r) D - 1 2w k 5 D - 1 (k - k') (A.2) 
with u;£ = k 2 + m 2 , and introduce the Lorentz invariant measure 

d D - x k /A , 

dk = tn Nn - — . A.3 
(2 7 r) I) - 1 2w k v ' 

If the intermediate A^-particle state consists of momenta ft, the N-body phase space is 
defined by 

(N \ N 

Using these conventions we have for the dimensions of the 2^2 scattering amplitude, 
[T(s,t)} =M A ~ D . 
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If we now restrict ( |A.1| ) to forward scattering, e.g. a = (3, we can replace the LHS by 
2% ImT(s, 0), and on the RHS we recognize the sum of the square of the amplitudes which 
enters in the definition of the total cross section. Recall that this is defined as 



a T = a (a — > all) 



1 



4a/ (pi ■ P2) 2 - m\m\ 



2n) D J2 I d$> N \T aN \ 2 . (A.5) 

N 



Notice that the prefactor in square brackets goes to l/^-Ei-E^) when s ^> m^m^. We 
are now ready to state the optical theorem, which is nothing but a direct consequence of 
unitarity: 



ImT(s, 0) = 2yJ {pi ■ p 2 ) 2 - m\ m z 2 a T (s) -> sa T (s) . (A.6) 

We can also relate the coefficients in the partial wave projections (|2.6|) , where the optical 
theorem takes the form (in the s ^> m 2 limit) 



Imf l (s) = 8(2ir) 2D - 2 ^y~ D/2 J2S D (PN-Pc<)\fi(s,{N})\ 2 , (A.7) 

N 



from which (|2.12| ) follows. In this expression the fi(s, {N}) are the partial wave projections 



of the generic intermediate states, considered modulo an overall rotation. The sum runs 
over all possible such subclasses of states [0. Performing the sum over / on both sides 
reproduces the optical theorem. 

As we emphasized in this paper, due to the masslessness of gravity we expect singu- 
larities at t = 0. We noticed before that the IR singularities can be removed by working 
in D > 4. From the definition of the cross section we promptly discover that we actually 
need even higher D for it to be well defined. This follows from the elastic cross section; 
( |2.5|) gives probability 

m 2 ~ p . (a.8) 

This Rutherford-like singularity is tamed for D > 6 by the integration over solid angle, 
with measure sin 0-3 9, giving a finite cross section. Once the cross sections are finite the 
optical theorem ( |A.6| ) shows that ImT(s, 0) is also finite. One may be tempted to push 
the partial wave expansion to t > 0, but this attempt fails once we realize that t = is 
indeed also a threshold for graviton production, and the partial wave expansion will not 
converge past that point. The finiteness of ImT(s, 0) is due to the fact that in higher 
dimensions the threshold behavior scales as a power of momentum, e.g. ~ (— t) a , rather 
than logarithmically as we are used to encountering in four dimensional field theories. This 
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is intimately linked to the softness of the IR divergences in D > 4 due to the promotion 

,4 ,D 

of the measure in the loop integrals from p^yr to / 27r yr> ' ^ * s then easy to see that the 
expansion of the derivatives of T(s, t) at t = will not converge and we cannot analytically 
continue the partial wave decomposition to positive values of t. 

A final comment is in order. The reader may be puzzled by the fact that the Born 
approximation in (275) seems to have a divergent imaginary part as t — > from the ie 
prescription. A careful analysis shows that is indeed not the case, and such singularity 
only arises in the plane-wave limit and disappears as soon as we take into account wave 
packets. The real part of the amplitude is large, but finite, and give rise to a finite 
contribution in the cross section as in (|A.8|) . 
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